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High-Mass Starless Cores 

T.K. Sridharan 1 ' 4 , H. Beuther 1 ' 4 , M. Saito 2 ' 4 , F. Wyrowski 3 - 4 , P. Schilke 3 ' 4 

ABSTRACT 

We report the identification of a sample of potential High-Mass Starless Cores 
(HMSCs). The cores were discovered by comparing images of the fields containing 
candidate High-Mass Protostellar Objects (HMPOs) at 1.2mm and mid-infrared 
(8.3/im; MIR) wavelengths. While the HMPOs are detected at both wavelengths, 
several cores emitting at 1.2mm in the same fields show absorption or no emission 
at the MIR wavelength. We argue that the absorption is caused by cold dust. 
The estimated masses of a few 10 2 M Q - 1O 3 M and the lack of IR emission 
suggests that they may be massive cold cores in a pre-stellar phase, which could 
presumably form massive stars eventually. Ammonia (1,1) and (2,2) observations 
of the cores indicate smaller velocity dispersions and lower rotation temperatures 
compared to HMPOs and UCHll regions suggesting a quiescent pre-stellar stage. 
We propose that these newly discovered cores are good candidates for the HMSC 
stage in high-mass star-formation. This sample of cores will allow us to study 
the high-mass star and cluster formation processes at the earliest evolutionary 
stages. 

Subject headings: stars: formation - stars: massive - ISM: dust - ISM: clouds 



1. Introduction 

The area of high-mass star-formation is witnessing significant progress. Systematic 
studies have uncovered several High-Mass Proto-Stellar objects (HMPOs) in a pre-Ultra 
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Compact Hn (UCHu) region phase (Molinari et al 1996,2002; Sridharan et al 2002; Beuther 
et al 2002a). The ubiquity of the outflows in these objects points to essential similarities 
between the high- and the low-mass star-formation processes viz., the presence of a disk- 
accretion phase (Zhang et al 2005; Beuther et al 2002b). However, the characteristics of the 
High-Mass Starless Core (HMSC) stage which must precede the HMPO stage has barely been 
studied (Evans et al 2002). In contrast, recent studies of low-mass starless cores have begun 
to constrain the initial conditions for low-mass star-formation (Motte et al 1998; Bacmann 
et al 2000; Alves, Lada & Lada, 2001; Harvey et al, 2003). Similar studies of HMSCs are 
important in answering one of the central questions in star-formation: how do star-formation 
processes produce the stellar Initial Mass Function? Since high-mass stars form in clusters, 
and a major fraction of all stars form in such environments, the structure of the HMSCs 
represents the initial conditions relevant for star-formation in general. 

In this letter, we present and study a sample of potential HMSCs identified by com- 
bining images at millimetre and mid-infrared wavelengths of HMPO fields. The candidate 
HMSCs belong to the general class of infrared-dark clouds (IRDCs), studied recently by 
others (Perault et al 1996; Egan et al 1998; Carey et al 2000, Bacmann et al 2000; Teyssier, 
Hennebelle & Perault 2002, Garay et al 2004). At the outset, we want to clarify two points: 
we use the terms starless and pre-stellar to indicate the lack of massive star-formation (> 8 
M Q ; spectral types B2 or earlier); the objects discussed here are candidates for the HMSC 
stage although this is not mentioned every time for brevity. 



2. Source Sample 

The data used to identify the HMSCs consist of continuum images at 1.2mm and 
8.3/um wavelengths, of a sample of 69 HMPOs we have been studying (Sridharan et al 
2002, Beuther et al 2002a,b, Williams, Fuller & Sridharan, 2004, 2005). The 1.2mm im- 
ages obtained with the MAMBO array at the IRAM 30M telescope have been described 
in detail by Beuther et al (2002a). The 8.3yum images are the archival data from the 
Mid-Course Space Experiment (MSX; Egan et al 1998), available at the MSX web site 
(http:/ /www. ipac.caltech.edu/ipac/msx/msx. html). 

The 1.2mm images, each covering a ~ 5'x 5' field around an HMPO, detected multiple 
sources in a number of fields with an average multiplicity rate of 2.2 (Beuther et al 2002). 
We argued in Sridharan et al (2002) and Beuther et al (2002a) that the objects in our HMPO 
sample are likely to be in various stages of evolution. This is expected, given the clustered 
mode in which high-mass star-formation is thought to occur. Further, it is conceivable that 
an earlier starless stage in the high-mass star-formation process should occur in the same 
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fields. With this motivation, we compared the 1.2mm and the MSX A band (8.3/im) images 
of the HMPO fields. Strong emission at the two wavelengths was found to be typical of the 
HMPOs, within reasonable position uncertainties. Significant position offsets between MSX 
and mm emission are accounted for in terms of a cluster setting and varying evolutionary 
stages (Sridharan et al 2002). However, in many fields, we also discovered 1.2mm emission 
features associated with either MSX absorption or no detectable emission. Figure 1 shows 
three examples. The strongest mm emission near the center of the field in Figure la is 
associated with the HMPO-IRAS18385-0512 which is also detected at 8.3/xm. The mm 
emission feature seen to the south east of the HMPO is associated with absorption at 8.3/rni. 
Figures lb (HMPO-IRAS18223-1243; also studied by Garay et al 2004) and lc present more 
examples, to indicate the range of morphologies, showing filamentary structures and a field 
of multiple cores where a stellar cluster may be forming. Following similar observations of 
low-mass starless cores (for example, Bacmann et al, 2000), we suggest that the absorption 
is caused by the same cold dust emitting at 1.2mm. The similarity of the morphologies of 
the emission and absorption supports this hypothesis. In the case of IRDC-18385-0512-3, 
shown in Figure la, using a temperature of 15 K (see section 3.1), and the 1.2mm flux, a 
mass of 344 M is estimated, making it a potential High-Mass Starless Core. 

We have identified a total of 56 candidate HMSCs, listed in Table 1 along with kinematic 
distances, 1.2mm fluxes and mass estimates. Each object is labelled with the HMPO field 
and the mm-emission clump number from Beuther et al (2002a). For completeness, we have 
included 26 prominent MSX absorption objects either lying outside the 1.2mm images or 
where the images are too noisy, extending the numbering. For mass estimates (30 objects), 
as for the HMPOs (Beuther et al 2002a, 2005; Hildebrand 1983), we used a value of 2 for 
the spectral index of dust emissivity (3, 100 for the gas to dust ratio, 0.1 /xm for grain size 
and 3 g cm~ 3 for grain mass density (k = 0.4 cm 2 /gm of dust), in addition to a temperature 
of 15K (section 3.1). In case of distance ambiguity, the near distance has been used - at the 
far distance MIR absorption is unlikely. As seen in the table, the cores are massive, with 
masses in the range of a few 1O 2 M - 1O 3 M . 

An examination of the High Resolution processed IRAS images (HIRES; Aumann, 
Fowler, Melnyk 1990) of the fields failed to detect any emission from dust heated by star- 
formation activity. Based on the Arcetri Catalog of water masers (Valdettaro et al 2001), 
the General Catalog of 6.7GHz methanol masers (Pestalozzi, M.R., Minier V. and Booth, 
R.S., 2005) and interferometric maser observations of the HMPO fields (Beuther et al 
2002c), we find two cases - IRAS18102-1800-1 & IRAS18151-1208-2 - to be associated with 
methanol/ water maser emission respectively, implying star-formation in these cores. This 
suggests that the cores presented here contain objects in multiple evolutionary stages. How- 
ever, all these stages represent phases before the pre-UCHll/HMPO phase, containing either 
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low-mass young stellar objects or early pre-cursors to the HMPOs. A further implication is 
that the maser activity may turn on very early in the star-formation process, and a majority 
of the cores presented here may be in a younger phase. 

3. Observations & Results 

To characterize the state of the dense gas in the sample, a subset of 34 objects from 
Table 1 was observed with the Emesberg 100-m telescope for NH 3 emission in Dec 2002. 
The (1,1) and (2,2) inversion lines were observed simultaneously in frequency switched mode 
with a velocity resolution of Av = 0.25 kms -1 and ~ 4 minutes of integration time per 
position. The system temperatures were mostly below 50 K. Both lines were detected in all 
the observed sources. For 18 sources with good quality spectra, rotation temperatures were 
derived following Ungerechts et al. (1986). For temperatures below 20 K, as here, kinetic 
temperatures are only marginally higher (Danby et al., 1988). Figure 2 shows an example 
NH 3 (1,1) spectrum, for the object IRAS18385-0512-3 (whose image is shown in Figure la). 
This core has a small linewidth of 1.1 kms -1 . 



3.1. Ammonia line-widths and temperatures 

Figure 3 compares the distributions of the ammonia line widths and rotation temper- 
atures for the HMSCs, HMPOs and UCHn regions (Sridharan et al 2002, Churchwell et 
al, 1990). The HMSCs exhibit a smaller line width with a median of 1.5 kms -1 (mean 1.6 
kms -1 ) compared to 1.9 kms -1 (mean 2.1kms~ 1 ) for the HMPOs. This suggests smaller in- 
ternal motions and therefore more quiescent conditions, supporting the pre-stellar nature of 
the HMSCs. The rotation temperatures indicate colder conditions in HMSCs with a median 
value of 16.9 K (mean 15.3) compared to 18.5 K (mean 22.5) for HMPOs, again suggesting 
more quiescent conditions in the HMSCs. The corresponding ammonia line widths and tem- 
peratures for UCHn regions are 3 kms" 1 and 22 K (Churchwell et al 1990). Thsese results 
support the picture where the objects evolve from HMSCs to HMPOs to UC-Hn regions. 

Although the ammonia rotation temperatures for the HMSCs are only marginally lower 
than for the HMPOs, where grey-body dust temperatures were estimated to be much higher 
(Sridharan et al, 2002), we believe that the dust in the HMSCs is cooler, and well character- 
ized by the gas temperature. For the HMPOs, the best fits to the spatial and the spectral 
energy distributions at FIR and mm- wavelengths imply temperatures increasing inwards 
(Williams, Fuller & Sridharan, 2005). Ammonia emission traces the cooler outer regions of 
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the cores, leading to lower rotation temperature estimates. In contrast, for externally heated 
cores, as suggested here for the HMSCs, the dust temperature is expected to be similar to 
the gas temperature in the outer regions or may even be marginally lower (Li, Goldsmith 
& Menten, 2003). The lack of FIR emission towards the HMSCs - in cases where the the 
dominant HMPO in the field is not nearby, no emission is seen even at 60fim or 100/im 
wavelengths in the HIRES images - indicates the absence of dust heated by an embedded 
object. 

The distributions of ammonia rotation temperatures and linewidths of the HMSCs and 
HMPOs were subjected to a Kolmogorov-Smirnov test, resulting in 0.90 and 0.62 for the 
probabilities that they are different. The masses of HMSCs are distributed over a similar 
range as the HMPOs. Virial masses estimated from ammonia line widths and 1.2mm sizes 
are lower, but comparable to dust derived masses, consistent with the HMSCs being in 
equilibrium. Because of the many uncertainties in mass estimates in both the cases, we do 
not make further comparisons. 

As already noted, we use the term starless here to imply the absence of a massive star. 
From the IRAS flux limits, a luminosity upper limit of ~100 L Q is derived for a distance 
of 4 kpc, corresponding to early A spectral types (~ 2-3 M ). Assuming a star- formation 
efficiency of 30%, an HMSC mass of 550 M (sample median), and an IMF with power-law 
indices 2.3 and 1.3 for masses > 0.5 M and 0.08 - 0.5 M respectively (Kroupa, 2004), we 
estimate that one star of mass > 20 M could form in an HMSC. While there is no current 
evidence for massive star-formation in these cores from their infrared fluxes, more sensitive 
seraches or studies like the Spitzer GLIMPSE survey will presumably uncover lower mass 
star- formation or hitherto unknown early stages of high-mass star-formation. 

In summary, we have identified a sample of candidate high-mass starless cores, likely to 
contain objects in multiple stages of evolution prior to the formation of massive stars. The 
sample provides an opportunity to determine initial conditions for high-mass star/cluster 
formation and star-formation in general. The study of their internal structure at multiple 
wavelengths using instruments like the Sub- millimeter Array (SMA) and the Spitzer telescope 
will pave the way for the understanding the origin of the stellar Initial Mass Function (IMF), 
a central puzzle in star-formation. 

H.B. thanks the Emmy-Noether-Program of the Deutsche Forschungsgemeinschaft for finan- 
cial support (DFG, grant BE2578/1). 
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Fig. 1. — MSX A-band (8/im) images (black is bright) with 1.2mm emission contours: (a) IRDC- 
18385-0512-3; (b) and (c) show fields 18223-1243 and 18454-0158 with filaments and multiple cores. 
The star symbols mark cores lacking good 1.2mm measurements. 
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Fig. 3.— The cumulative distributions of NH 3 (1,1) line widhts for HMSCs and HMPOs with the 
median for UCHii regions shown as a vertical line (left). The cumulative distributions of NH3 
rotation temperatures for HMSCs and HMPOs (right). 
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Table 1. Table 1 - Candidate High-Mass Starless Cores 
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Table 1 — Continued 
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56 


4.7 d 


0.4 


51.9/99.6 


1.5/1.9 




549 


41 


IRDC-IRAS18454-0158-3 


18 47 55.8 


-01 53 


34 


6.0 a 


0.7 


93.7/97.6 






1565 


42 


TRDC-TRAS18454— 01 ^8-5 


18 47 58.1 


-01 56 


1 


fi o a 


3 


Q3 9 


1 5 


1 7 fi 


fi71 


43 


IRDC-IRAS18454-0158-8 


18 47 50.5 


-01 57 


28 






94.6 








44 


IRDC-IRAS18454-0158-9 


18 47 50.0 


-01 53 


46 






95.6/101.8 








45 


IRDC-IRAS18454-0158-10 


18 48 10.0 


-01 53 


29 






100.4 








46 


IRDC-IRAS18454-0158-11 


18 48 07.0 


-01 53 


25 














47 


IRDC-IRAS18454-0158-12 


18 48 05.7 


-01 53 


28 














48 


IRDC-IRAS18460-0307-3 


18 48 36.0 


-03 03 


49 


5.2 


0.3 








504 


49 


IRDC-IRAS18460-0307-4 


18 48 46.0 


-03 04 


05 


5.2 


0.3 








504 


50 


IRDC-IRAS18460-0307-5 


18 48 47.0 


-03 01 


29 


5.2 


0.5 








840 


51 


IRDC-IRAS18530+0215-2 


18 55 29.0 


+02 17 


43 














52 


IRDC-IRAS19175+1357-3 


19 19 52.1 


+14 01 


52 


l.l a 


0.2 


7.7 


1.5 


15.6 


15 


53 


IRDC-IRAS19175+1357-4 


19 19 50.6 


+14 01 


22 


l.l a 


0.2 


7.7 


1.5 


13.2 


15 


54 


IRDC-IRAS19410+2336-2 


19 43 10.2 


+23 45 


04 


2.1 


1.4 


21.4 


2.0 


18.3 


383 


55 


IRDC-IRAS20081+2720-1 


20 10 13.0 


+27 28 


18 


0.7 


0.6 


5.5 






18 


56 


IRDC-IRAS22570+5912-3 


22 58 55.1 


+59 28 


33 


5.1 


0.5 


-47.1 


2.6 


16.2 


807 



a different velocities, kinematic distances and masses compared to those in Sridharan ct al (2002) and Buether 
ct al (2002) 

b a dust temperature of 32.7 K was used to compute this mass; significantly higher temperature than the rest 
of the sample 

c thc RA reported in Bcuther ct al (2002) for this object is in error by ~30" 
d an average of the kinematic distances for the two velocity components 
c may not qualify as high-mass; included for completeness 



